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The nuclear density of states plays an important role in nuclear reactions. At high energies, above 
a few MeV, the nuclear density of states is well described by a formula that depends on the smooth 
single particle density of states at the Fermi surface, the nuclear shell correction and the pairing 
energy. In this paper we present an analysis of the low energy behaviour of the nuclear density 
of states using the saddle point approximation and extensions to it. Furthermore, we prescribe a 
simple parabolic form for excitation energy, in the low energy limit, which may facilitate an easy 
computation of level densities. 
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One of the important ingredients in the Hauser-Feshbach approach to the calculation of nuclear reaction rates 
important for astrophysical interest is the nuclear density of states |lj] . In fact uncertainties in the nuclear density of 
~Z?\ ' states is a leading cause of errors M in these calculations. 

The studies of nuclear level densities dates back to the 1950s, with work by Rozenweig M, and Gilbert and 
Cameron pLEJ . The usual technique is to calculate the partition function and then invert the Laplace transform using 
the saddle point approximation. At energies sufficiently high for shell and pairing effects to be washed out the density 
of states is given in terms of the single particle density of states at the Fermi surface (and its derivatives) the shell 
correction energy and pairing energy. Most statistical model calculations use the back shifted fermi gas description [0 . 
Monte-Carlo shell model calculations || as well as combinatorial approaches || show excellent agreement with this 
phcnomcnological approach. At lower energies the results are more problematic and typically crude extrapolations 
t-H . from the higher energy are used. 

In this paper we study the nuclear level density with an emphasis on the lower energy region, using a single particle 
shell model. The dependences in the two regimes are rather different. In contrast to higher energies, where the density 
of states depends on the shell correction and the smooth single particle density of states, in the lower energy regime 
the density of states depends on the separation of single particle levels and their degeneracy. Moreover, at very low 
energy the saddle point approximation itself breaks down. We show that in this region the correction suggested by 
Grossjean and Feldmeier [J7 gives dramatic improvements. 

In the next section we review the saddle point approximation in the context of nuclear level density calculation. We 
use thermodynamic identities to rewrite the equations in simpler form compared to the usual ones 0. Furthermore, 
we discuss the possible ways to simplify the evaluation of level densities at low energies. A temperature dependent 
parabolic equation for excitation energy seems to be a good choice. The corrections suggested m and the corresponding 
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modifications to the equation are discussed in section 3. Finally in section 4 we discuss our calculation and the results. 



II. THE SADDLE POINT APPROXIMATION 



The grand canonical partition function for two type of particles can be written as: 

e n = Y^ exp(a N N' + a z Z' -f3E') (2.1) 



N',Z',E' 

where the sum is over all nuclei with TV' neutrons, Z' protons and over all energy eigenstates E 1 . r = (3~ l , is the 
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temperature and Hn(z) — ° w i z) is the chemical potential for neutrons (protons). The sum over eigenstates can be 



substituted by an integral: 
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e n = 



= Y^ I dE' P (E',N',Z')cxp(a N N' + a z Z' - [3E') 
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(2.2) 



where p(E',N',Z') is the nuclear density of states. It represents the density of energy eigenvalues for the nucleus 
(TV', Z') at the energy E'. The above equation also shows that the grand partition function can be considered a 
Laplace transform of the nuclear density of states. The inversion integral is: 
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where the entropy S = fl — ajqN — azZ-\-(3E . The above contour integrals are also known as Darwin- Fowler integral. 
This integral is usually done by the saddle point approximation and in this section we will explore this approximation. 
The location of the saddle point is defined by the equations, 
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The path of the integration can be chosen to pass through this point. By expanding the exponent S in Taylor 
series about the saddle point and retaining only the quadratic terms, the nuclear density of states in the saddle point 
approximation can be written as: 
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where D is the determinant of the second derivative of S with respect to the parameters oin , olz , and f3. The 
determinant can be simplified to a product of factors by changing the variables which are held fixed when the 
derivatives are performed. 

The determinant is written as: 
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To simplify the above determinant we change the independent variables to r = 1//3 , /j,n = raw , and /xjv = tolz 
and change the dependent variable to ft' = tVL = tS + [InN + P zZ — E. In terms of the new variables the equations 
determining the saddle point are: 
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Using these equations the determinant can be written as: 
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In deriving this result we have used the fact that in a determinant addition of a multiple of row (column) to another 
row (column) does not change the value of the determinant. In the first row the derivatives are at constant Hn and 
[iz ; in the second at constant r and /iz, and the third at constant r and /xjv • The variables which are held constant 
can be changed using the following equations: 
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and 
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Subtracting 4% I „ times the second column and 4# I „ times the third column from the first column we have: 
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This procedure can be repeated to yield: 
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Note the progression on the variables held constant. This procedure can be extended in an iterative manner to any 
number of constants of motions. The main gain is that we now have a simple product rather than a determinant. 

Let us now try to understand how this new form may help us in practice. To calculate the density of states as 
a function of energy for fixed particle number we need the entropy S as a function of energy at hxed N and Z for 
the exponent in the numerator. The temperature can then be obtained from dS/dE\ NZ — 1/r and dS/dr\ NZ = 
— l/(T 3 d 2 S/dE 2 ). This leaves the derivatives of the particle numbers to be separately evaluated. Thus we have three 
independent functions to parameterize. 

To see that this modified form of the density of states agrees with the standard form we consider the in dep endent 
particle model where with a consta nts si ngle particle density of states g. The entropy then is S = 2\faE where 
a = ng/6. The temperature is r = y/E/a, dS/dr\ NZ — a/2, and dN/dfiN — dN/d/xz — 9- This then gives the well 
know formulae: 
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Now let us consider a normal quantum system with a discrete spectrum. In this case there are in general no closed 
forms for the various functions, so we consider the zero temperature limit of dN/dfijy. We start with an open shell 
situation where there is a partially filled shell. For this discussion it is only necessary to consider the properties of the 
partially filled level. We take the level to have a degeneracy of g\ , an energy of t\ and a partial occupancy of d. As 
t goes to zero the saddle point condition for the number of particle becomes <?i/(l + exp[(ei — li)/t]) = dg\. Solving 
for fi we have fi = e\ + rln[d/(l — d)]. Note that this formulae breaks down for d equals zero or one corresponding to 
closed shells. The derivative dN/dfx is given as g\d(l — d)/r. Note that it diverges as r goes to zero. 

For a closed shell it necessary to consider two levels, the last filled level and the first unfilled level. We denote 
the energies and degeneracies of these levels as ei,gi and £2,32- The saddle point condition is now 32/(1 + exp[(ei — 
A 4 )/ 7 "]) + 32/(1 + exp[(e2 — m)/ t D = <7i- Solving for [i we have 11 = (ei + C2)/2 — T\a[g\/g2\ for small temperatures. 
The derivative dN/dLi is given as 2^/gig 2 exp[(ei — e2)/(2r)]/r. This goes to zero exponentially fast as r goes to zero. 
Note that in neither of the cases above shell correction is involved. 

The above discussion is useful since S is a function of the energy. Here again one may use a few trick. It turns out 
to be easier to parameterize E as function of r. Since rdS = dE, one can write S — J„ -p-jp-dr' + S(t = 0). The 
last term is the integration constant and is given once the degeneracy of the ground state is known. It contributes to 
exponent but not to the denominator where derivatives are taken. 

Next one needs £ as a function of r. For many systems there is a quite reliable approximation. For very low 
temperature, much less then the level spacing, the energy does not change significantly. However above some critical 
temperature ,r , it starts to increase rapidly. For temperature nears this region the energy can be parametrized quite 
simply by 
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We have checked this approximation using a simple shell model and found that it works quite well except if there 
are more then one level approximately equal distant from the Fermi surface. The parameters To and c depend on the 
level spacing and degeneracy near the Fermi surface. Again they do not depend on the shell correction. 

Before being useful at very low energies a short-coming of the saddle point approximation must be overcome. It is 
well known that at low energies the saddle point approximation tends to diverge as the denominator goes to zero, In 
many cases this problem can be fixed by using a technique from ref. M which handles the contribution to the nuclear 
density of states from the ground state delta function explicitly. 



III. MODIFIED SADDLE POINT 

In ref. H Grossjean and Feldmeier have proposed a modification of the saddle point method to remove the diver- 
gences of the level density at the ground state. Introducing explicitly the ground state energy E g (A) as the lower 
boundary, the density of state becomes 



ui(E*,A) = lu{E* + E g (A), A) - 5(E*)S{A - A ) 



(3.1) 



where E* = E — E g is the excitation energy and Aq is the mean particle number. The corresponding modified grand 
canonical potential is given by, 
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The chemical potentials for neutrons and protons are given by nz = ^g- and [In = ^ff- respectively, do being the 
ground state occupancy. 

The nuclear level density in the modified saddle point approximation becomes, 
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where D is the determinant in the form of eq. [2.8| , defined in terms of S, the entropy corresponding to the potential 
fl. The steps of the previous section (eq.(|2.8|) - eq.(2.13|)) can be retraced and then one gets, 



D = -r b 



dS_ 



dN 

NZ 


dZ 
rz d^ z 



(3.5) 



The derivation of eq.(3.5) is straightforward as it depends only on the thermodynamic relations of the quantities 
involved and not on their explicit forms. 



The computation of the level density using eq.(3.4) will depend on the relations between the usual and the modified 
thermodynamic quantities as the usual quantities are directly related to the single particle shell model states. 
The modified saddle point conditions in terms of the usual thermodynamic potential becomes, 
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where E* = E*(l — Y) (see eq.( |3.3| )). For a = ajv(z)i A — N(Z). The derivatives of the entropy and numbers are 
related as, 
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Using the above equations one can compute the level density (eq.(3.4)) in terms of single particle states. 



IV. DISCUSSION 

The single particle energies, required for the evaluation of different thermodynamic quantities are obtained for the 
Nilsson shell model. The values of the constants associated with l 2 and l.s are taken from ref. ||. Here it should be 
mentioned that the level densities are strongly dependent on the single particle energy levels. Hence for astrophysical 
applications one should make a judicious choice for the model as well as the constants. 

We first calculate the level densities for different nuclei, for both the usual as well as modified saddle point approx- 
imations, using all the filled and a equal number of unfilled levels. It should be noted the for low excitation energies 
(of the order of first excitation level or less) the inclusion of only the last filled and the first filled level, as described 
in section 2, is sufficient for the evaluation of the level densities. 



A comparison of the level densities from cq.(p.3|) and eq^^) for nuclei 32 S, S8 Sr and 20S Pb are shown in fig. 1, 
fig. 2 and fig. 3 respectively. The modified saddle point results are shown by curve (a) and the usual saddle points 
results are shown by curve (c) in the above figures. As evident from the graphs the usual saddle point does show a 
divergence at the low energies whereas modified version goes to zero smoothly. This is due to the fact that entropy in 
the modified version goes to zero much faster as it takes into account the nonavailability of states below the ground 
state energy. Moreover, at low energies the differences are more pronounced for the lighter nuclei like 32 5* compared 
to 208 Pb. The differences can be attributed to the fact that in modified prescription the thermodynamic potential 



gets an additive contribution compared to the usual one as shown by eq.(3.2). 

As discussed in section 2. we try to fit the excitation energies with parabolic form as given in eq.(2.14). These fits 
for different nuclei are shown in fig. 4, fig. 5 and fig. 6 respectively. In figures (4-6) we have plotted the variation of 
excitation energy with temperature. For low energies, the fitted value is in good agreement with exact values. Next 
we calculate the entropy and its derivative, using the steps given in section 2, from this fitted expression for excitation 
energy, the derivatives of N and Z being the same as in preceding paragraph. Using these we calculate the level 
densities for different nuclei. A comparison of the level densities from the full calculation from modified saddle point 
approximation (curve (a)) and the one using fitted excitation energies (curve (b)) are shown in fig. 1, fig. 2 and fig. 
3. It is obvious from the graphs that the fitted excitation energy gives a better agreement with exact calculation for 
heavier nuclei. 

To conclude, we have shown that the modification of the saddle approximation is necessary for the correct evaluation 
of the level densities at lower energies. One can simplify the equations substantially using the thermodynamic 
identities. Furthermore, a parabolic prescription for the excitation energy may be useful for easier computation of 
the level densities. More work in this direction is needed to make the methodology useful for direct application to 
astrophysical reactions. 
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FIG. 1. Level density for S; (a) modified saddle point, (b) corresponds to the fitted excitation energy as in fig. 4 and (c) 
usual saddle point. The vertical dashed line gives the position of 1st excitation level 
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FIG. 2. Level density for Sr; (a) modified saddle point, (b) corresponds to the fitted excitation energy as in fig. 4 and (c) 
usual saddle point. The vertical dashed line gives the position of 1st excitation level 
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FIG. 3. Level density for 208 Pb; (a) modified saddle point, (b) corresponds to the fitted excitation energy as in fig. 4 and 
(c) usual saddle point. The vertical dashed line gives the position of 1st excitation level 
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FIG. 4. Excitation energy for 32 S; (a) actual excitation energy nd (b) corresponds to the fitted form eq.(2.14) with £o=5.62 
and t =0.11 
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FIG. 5. Excitation energy for Sr; (a) actual excitation energy nd (b) corresponds to the fitted form eq.(2.14) with Eq = \§.71 
and r =0.14 
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FIG. 6. Excitation energy for 208 Pb; (a) actual excitation energy nd (b) corresponds to the fitted form eq.(2.14) with _Eo=40.0 
and ro=0.20 
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